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Abstract 

We calculate fermionic response in domain wall backgrounds of four-dimensional gauged 
supergravity interpolating between distinct stable AdS vacua. The backgrounds, found by 
Bobev et ah, are holographically dual to zero-temperature states of ABJM theory at finite 
density for monopole charge and are similar to zero-temperature limits of holographic su¬ 
perconductors, but with a symmetry-breaking source as well. The condensed scalar mixes 
charged and neutral fields dual to composite fermionic operators in the Dirac equations. 
Both gapped and gapless bands of stable quasiparticles are found. 
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1 Introduction 


Understanding the behavior of strongly coupled fermionic systems at nonzero density is of 
great interest. Holographic methods using the AdS/CFT correspondence are a powerful tool 
for exploring such systems, and they have been employed prohtably both from a “bottom-up” 
perspective, where a custom gravity theory can be tailored to produce the desired dynamics, 
and a “top-down” perspective, where helds and solutions of supergravity and string theory 
can be matched precisely to operators and states of known dual quantum held theories. 
Nonzero density means turning on a chemical potential for a conserved charge, associated 
to a background gauge held on the gravity side. 

One class of solutions that has been studied extensively leaves the corresponding symme¬ 
try unbroken. Holographic Fermi surfaces at zero temperature were studied, initially from 
a bottom-up perspective of generic fermions in Reissner-Nordstrom backgrounds [DEI Em 
and later from a top-down perspective in gravity duals to four-dimensional A/” = 4 Super- 
Yang-Mills theory and three-dimensional ABJM theory involving more complicated black 
hole geometries with running neutral scalars lainiiDEiEiEniin]. The game is to calculate 
retarded fermionic Green’s functions by solving the Dirac equation in the appropriate gravity 
background, interpreting bulk fermion (quasi)normal modes at zero energy as holographic 
Fermi surfaces and studying the corresponding dispersion relations around them. In bottom- 
up models, parameters of the bulk Lagrangians can be adjusted so that the excitations near 
the holographic Fermi surface resemble those of either a Fermi liquid (with stable quasi¬ 
particles) or a non-Fermi liquid (without stable quasiparticles). In the top-down models 
of strongly coupled A/" = 4 Super-Yang-Mills and ABJM theories, however, the bulk La- 
grangian is hxed, and the fermionic excitations that appear reflect actual dynamics of these 
held theories at hnite density, modulo possible instabilities. Generic backgrounds in both 
cases have nonzero entropy at zero temperature, and associated dispersion relations char¬ 
acteristic of non-Fermi liquids. At special values of the chemical potentials the zero-point 
entropy vanishes and an energy band of absolutely stable quasiparticles appears around the 
Fermi surface. Fermions at hnite temperature were also studied for top-down models, leading 
to agreement with the zero-temperature cases. 

Another interesting class of backgrounds, which we focus on here, breaks the symmetry 
of the conserved charge. This includes the holographic superconductors, in which a charged 
scalar condenses outside the horizon in the gravity background for sufficiently low temper¬ 
ature HD ESI m US]. The zero-temperature limits of such backgrounds are expected to 
be horizonless domain wall-type geometries, where Lorentz invariance or even conformal in- 
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variance may be regained in the infrared dSl HZl HE]. Snch backgronnds conld be related 
to the non-Fermi liqnids discnssed previonsly, as real-world non-Fermi strange metals in 
high-temperatnre snpercondnctors become hidden behind a snpercondncting dome below a 
critical temperatnre, bnt the zero-temperatnre qnantnm critical point is still thonght to con¬ 
trol their behavior. Symmetry-breaking systems have been stndied nsing bottom-np generic 
fermion actions in [19]. In 121 . a Majorana Ynkawa conpling was stndied where a charged 
fermion conples to itself (not its conjngate) as well as a “Cooper pair” scalar with twice the 
charge. It is interesting to ask whether the Fermi system becomes gapped in the presence of 
symmetry breaking; [20] fonnd that the Majorana conpling conld lead to the generation of 
snch a gap by a “level crossing” repnlsion between two lines of poles in the Green’s fnnction. 

Natnrally, it is interesting to consider symmetry-breaking systems from a top-down per¬ 
spective. Holographic snpercondnctors descending from string/M-theory were constrncted 
in [211 [22112211211120], and one may ask abont the fermionic response of these systems, ns¬ 
ing the Dirac eqnations derived from the top-down theory. Here we will focns on a class 
of symmetry-breaking geometries described in [26]. These backgronnds are solntions to the 
eqnations of fonr-dimensional M = 8 ganged snpergravity invariant nnder a SO(3) x SO{3) 
snbgronp of the SO{8) gange gronp, and besides the metric involve a gange held correspond¬ 
ing to the chemical potential and a charged scalar. These geometries are domain walls, 
interpolating between the maximally symmetric AdS^ vacnnm in the UV and a different 
AdS^ region corresponding to a nonsnpersymmetric critical point of the scalar potential in 
the IR, and as snch are similar to the zero-temperatnre limit of holographic snpercondnc¬ 
tors. Both the UV and IR hxed points are known to be stable nnder pertnrbative scalar 
hnctnations 123- The dnal descriptions of these solntions are states of three-dimensional 
snperconformal ABJM theory with a chemical potential for monopole charge, and both a 
sonrce and an expectation valne tnrned on for charged operators dnal to the bnlk scalar; as 
a resnlt the system is not precisely a snpercondnctor, since the associated U{1) is broken 
explicitly as well as spontaneonsly. 

To nnderstand fermion response in these ABJM states, we analyze the spectrnm of 
fermionic hnctnations in the holographically dnal backgronnds. We stndy spin-1/2 modes 
whose SO{3) x SO{3) qnantnm numbers prevent them from mixing with the gravitino helds. 
While they do not mix with the gravitino, these spin-1/2 modes mix with each other; another 
system of mixed fermionic excitations in a symmetry-breaking background was studied in 
[22]. We study a simplihed, non-chiral version of the full top-down fermion mixing matrix, 
where in a more intricate generalization of the Majorana “Cooper pair” coupling of 121 , 
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a charged fermion mixes with a neutral fermion via the condensed charged scalar, which 
in turn has a Majorana-type coupling to its own conjugate. Thus the analog of a Cooper 
pair in this system is a condensation of a charged/neutral bound state. Since the charge is 
monopole charge, the associated charged operators may be thought of as bound states of 
fundamental fermions with vortices, so-called composite fermions. 

We study two backgrounds, one having a source for a fermion bilinear and an expectation 
value for a boson bilinear, and the other with the roles reversed. We hrst elaborate on the 
analysis of [26] of the conductivities of these backgrounds. In geometries approaching AdS in 
the infrared, the emergent Lorentz invariance generates a light-cone structure on the energy- 
momentum, and modes living outside this light-cone have regular (rather than infalling) IR 
boundary conditions; normal modes in this region are associated to fluctuations in the held 
theory with zero dispersion. We determine the locations of such normal modes, and hnd in 
both cases two lines of modes, one gapped and the other ungapped. By studying the pole 
structure of the matrix of Green’s functions, we can see that each line is a mixture of both 
charged and neutral fermionic excitations. This matrix also provides information about the 
presence of unstable huctuations inside the lightcone. 

Following the ideas of [20], it is natural to inquire how the charged fermion/neutral 
fermion coupling in our system of Dirac equations affects the dynamics, and in particu¬ 
lar whether it causes a repulsion between lines of poles. We study a modihcation of the 
Dirac equations removing the couplings between different fermions, and see that without the 
charged/neutral coupling there is an ungapped, purely charged band, a purely neutral band 
that asymptotes to the origin in energy/momentum space, and a new gapped band with 
the conjugate charge; each pair of bands has a point of intersection. The charged/neutral 
coupling thus has a number of effects: it pushes the third band outside the region of stable 
excitations, it repels the crossing between the other two bands while mixing the charged and 
neutral contributions, and it turns what was the neutral band away from the origin leaving 
it fully gapped. While modifying the Dirac equations departs from the top-down structure 
of AT = 8 supergravity, it allows us to see how the couplings are responsible for the structure 
of the quasiparticle excitations. The fully top-down, chiral fermion mixing matrix is studied 
in [2E]. 
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2 Gauged supergravity and ABJM theory 


In this section we review the maximally snpersymmetric gauged supergravity theory in four 
dimensions and a particular truncation of it in which the backgrounds we will consider can 
be constructed, as well as its duality to (2+l)-dimensional ABJM theory. 


2.1 4D J\f = 8 gauged supergravity 

The four-dimensional maximally supersymmetric gauged supergravity theory [2^ [30] is the 
consistent truncation of eleven-dimensional supergravity compactified on a seven-sphere to 
retain only the supermultiplet of the four-dimensional graviton. The bosonic degrees of 
freedom are the vierbein, 28 gauge fields in the adjoint of the gauge group SO{S), and 70 
real scalars; the fermions are 8 Majorana gravitini and 56 Majorana spinors. The scalars 
parametrize the coset Ejij)/SU{^) as a 56-bein, which can be written 


V 


u. 


ij 


MU 


VijKL 


U 


kl 


KL 


exp 


0 (j)l.JKL 
^IJKL 0 


where the complex = 4>*ijkl obey the self-duality relation 


( 1 ) 


(pIJKL — -^^IJKLMNPQ ■ (2) 

In the second equality of ([^ we have gauge-fixed the internal SU{8) symmetry. This “unitary 
gauge” removes the distinction between 50(8) index pairs [IJ] and SU{S) pairs [ij], and 
allows us to associate definite 50(8) representations to all the helds: the scalars split into a 
35v of parity even scalars and a 35c of parity odd pseudoscalars, the gravitini are in the 8s, 
and the Majorana spinors are in the 56s. 

For the backgrounds we study we will be interested in a particular truncation of the 
gauged supergravity, retaining only modes invariant under an 50(3) x 50(3) subgroup of 
50(8). We will choose I = 3,4,5 and I = 6, 7, 8 to be the directions in which the two 50(3) 
groups act on the 8s; the 1,2 directions correspond to an additional 50(2) gauge symmetry 
commuting with 50(3) x 50(3). The truncation corresponds to an A/" = 2 gravity multiplet 
plus a hypermultiplet, with the bosonic sector consisting of the vierbein, a single graviphoton 
gauge held for the 50(2) in the 12-directions, and two complex scalars charged under it. 
Moreover, it is consistent with the equations of motion to set one complex scalar to zero, and 
we will do this in what follows. This 50(3) x 50(3)-invariant truncation is characterized 
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by the ansatz [HT] 


4>ijkl{x) = 


2^2 


cos a{x) + i Yjj^A sin a{x) i Zjj^^ 


(3) 


where A and a are four-dimensional scalars, and and are self dual (-I-) and anti-self 
dual (—) invariant four-forms on the scalar manifold, dehned as 


= 4! = 4! (4) 

ZJjkl = 4! {5]fL - SWL) = 4! . (5) 


The more general case of two complex scalars would involve four independent coefficients for 
the four tensors. 

Given the scalar ansatz ([^, it is a straightforward matter to obtain the Lagrangian of 
the truncated theory. The bosonic sector of the gauged JY = 8 theory in four dimensions 
can be written [30] 

= R- - 1 [F+„ - S‘A) FS" + h.c.] - 2P. (6) 

Here we have not yet imposed unitary gauge, and SO{8) and SU{8) indices are distinct. Let 
us discuss the terms in turn. The curvature scalar R is the usual Einstein-Hilbert term. The 
tensor Aijki determining the scalar kinetic terms follows from the dehnition 


D^V ■ V-' 


1 ( 0 \ 

Y Afxmnpq 0 j 


(7) 


where the covariant derivative acts on S'0(8) indices with the S'0(8) gauge held Aj/ and on 
SU{8) indices with the composite connection 


j~\ IJ o IJ I yj h IJ 


2^ , 


( 8 ) 


and a similar expression for . This expression also implicitly hxes the composite SU{8) 
connection B, which we will discuss in section The gauge helds have non-abelian held 
strengths of the standard form, with F+ = (F+ -|- iF^)/2 the 

self-dual part of the held strength; these couple to the scalars in their kinetic terms via the 
S'-tensor dehned as 

(««„ + ««") . ( 9 ) 
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The SU{8) covariant T-tensor 





jm 


KI 




( 10 ) 


in turn defines the tensors Ai and A 2 , 

jtij _ rpikj Ajkl _ _^nn[jkl] 

■^1 21 ^^ ^ 2 z 3 ^^ ^ \ 

which appear in the scalar potential, 

Evaluating the Lagrangian ([^ in the 50(3) x 50(3) invariant truncation gives 

e-^C = ^R- - d^Xd^X - _ gA^) {d^^a - gA^ - V, 

where has been set to one, the gauge field is A = A 12 , and the potential is 

V = ^{s^- Ss^ - 12 ) with s = sinh A. (14) 

It is easy to see that this potential has critical points at s = 0 and s = ±2. In anticipation 
of the domain wall geometry, we will refer to the corresponding values of the scalar A as 


( 12 ) 


(13) 


Auv = 0 


and 


Air = ± log(2 + 1 / 5 ), 


(15) 


corresponding to AdS 4 solutions with AdS radii Luv = and Lir = respectively. 


Solutions to the equations of motion coming from (13) will provide the classical backgrounds 


we wish to probe, and we will discuss them in more detail in section 


2.2 The holographic dual ABJM theory 

The maximally superconformal theory in three dimensions living on a stack of N coincident 
M2-branes is holographically dual to M-theory compactified on AdS^ x 5^; in the large- 
N limit this reduces to eleven-dimensional supergravity, and hence the four-dimensional 
gauged supergravity theory we have discussed describes a set of low-dimension operators in 
this theory. For a single M2-brane the theory is 8 free scalars in the 8v and 8 free spinors 
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in the 8c of the SO{8) R-symmetry; for > 1, however, the theory becomes interacting. 
While it can be characterized as the IR limit of three-dimensional Super-Yang-Mills theory, 
it is most conveniently formulated as ABJM theory. 

ABJM theory ( [32l |33l IMl ESI [36] ; for reviews see [371|38|) is a 3D U{N) x U{N) Chern- 
Simons theory at levels {k, —k) coupled to bifundamental matter. The manifest supersym¬ 
metry is W = 6 and the manifest global symmetry is SU{4) x U{l)b- For general k this 
represents the theory of N M2-branes at a Z*, orbifold singularity. However, for the cases 
fc = 1,2 there is an enhancement to W = 8 supersymmetry and SO {8) R-symmetry; the 
decomposition of the eight-dimensional representations of SO(8) into SU{A) x U{l)b are 

8v —>■ 4i © 4_i, 8c —>■ 4i © 4_i, 8s —6o © I 2 © 1 —2 • (16) 


Here we are interested in the k = 1 case, corresponding a stack of N M2-branes with no 
orbifold. 

The bifundamental matter may be written as four complex scalars A = 1... 4 in the 4 
of SU (4) and four complex spinors if^A in the 4; both sets of helds are in the N x N of f/ {N) x 
U{N) and neutral under U{l)b. Alone these helds do not assemble into complete 50(8) 
representations; however, they combine with monopole operators, representing the scalars 
dual to the gauge helds, into gauge-invariant objects with proper 50(8) transformation 
properties. We will denote by e'^’’ the monopole operator with 0(1);, charge q in the g-fold 
tensor product of N x N; monopole operators are neutral under 50(4). We then have 
gauge-invariant operators such as 

4i: Y^e^, 4_i: , 4i: , 4_i: , (17) 


assembling into complete 8v and 8c representations according to (16). It is these combi¬ 
nations that are analogous to the free bosons and fermions in the N = 1 case; the ABJM 
presentation fractionalizes the symmetry carriers into ordinary matter charged under 50(4) 
and monopole operators charged under 0(1);,, which bind into gauge-invariant “composite” 
bosons and fermions. 

The supergravity modes discussed in the previous subsection are dual to such gauge- 
invariant operators. These are described in the table, with the dual ABJM operators indi¬ 
cated schematically: 
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SUGRA Mode 




XlJK 

Re (j)ijKL 

Im cI)ijkl 

50(8) Rep 

1 

8 s 

28 

56s 

35v 

35c 

Dual ABJM Operator 


SI^ 


Y^/J 

y 2 

ijj'^ 

Gonformal dimension A 

3 

5/2 

2 

3/2 

1 

2 


The first three sets of operators are the energy-momentum tensor, supercurrents and 50(8) 
R-symmetry currents. The 28 R-symmetry current operators include 15 5t/(4) currents, one 
U{l)b current, and 12 additional operators including monopoles, corresponding to the 
decomposition 


28 —>■ 15o © Iq © 62 © 6—2 • 


(18) 


We note that while some operators include monopoles and some do not, the enhancement 
to full 50(8) symmetry means that they are all treated on equal footing. Indeed, one can 
imagine distinct embeddings of 5f/(4) x U{l)b inside 50(8) where a monopole operator in 
one case becomes a non-monopole operator in the other. 

Let us now discuss how the 50(3) x 50(3) x 50(2) subgroup of the previous section 
relates to the ABJM picture. Since one has a full 50(8) to work with, one can imagine 
embedding 50(3) x 50(3) x 50(2) in 50(8) in a way that does not play nicely with 
50(4) X 0(1);,. However, it is natural and convenient to take the simplest choice, where 
50(3) X 50(3) is realized as the subgroup of 50(4) = 50(6) under which 


6^(3,1)©(1,3), 4,4^(2,2), 


(19) 


and the remaining 50(2) of the gauge field (13) is simply 0(1)^ itself. Hence the charge 
carried by supergravity fields that condense in the backgrounds we will study next is most 
simply realized on the held theory side as the monopole charge. 


3 The Domain Wall Background 


The results of the previous section isolate a sector of maximal gauged SUGRA in O = 4 
whose bosonic content includes the vierbein, one 0(1) gauge held, and one complex scalar. 
The dynamics of this sector are encoded in the Lagrangian (13), and any solution to the 


corresponding equations of motion can in principle be uplifted to a solution of SUGRA in 
O = 11. This feature is particularly notable in that it will allow us to exploit the holographic 









duality between solutions of M-theory and states in ABJM theory to study the properties of 
an explicitly known held theory at strong coupling. 

We will be interested in solutions to the equations of motion wherein the scalar interpo¬ 


lates between the two hxed points (15) of P, and the geometry asymptotes to different AdS 4 
regions in the IR and UV. Holographically, these solutions are dual to zero-temperature 
states of ABJM theory in which both the low energy and high energy physics is governed by 
(distinct) conformal held theories. In |26j such domain wall solutions were found numerically, 
and their identihcation with zero-temperature limits of novel states in the ABJM theory was 
discussed in some detail. We now review these solutions, and provide further commentary 
on their dual holographic description. 


3.1 Interpolating gravity solutions 

The solutions of interest [26] are given in the ansatz. 


ds^ = —G{r)e -|- r^dx^ -|- 


dr^ 

wr 


A = J/(r)dt, A = A(r), a = 0. (20) 


This ansatz sets g = 1, equivalent to Luv = 1/v^- The equations of motion are 

,sinh2(2A)T2 

-r — 2rX — X , 

x/2 

+ ^ + + 


0 = — e^- 


2^2 

1 , P , ^^sinh2(2A)T2 ^ ^ 


0 ~ ^ A J- 

r2 G 4G2 


2G ’ 


0 = - 
0 = 


sinh^(2A)T 2T' 1 

sinh2(4A)T2 V , ,, [2 , G' 1 , ^ , ,, 
- U+G“2* 


( 21 ) 

( 22 ) 

(23) 

(24) 


Near the boundary, the solution approaches the maximally supersymmetric A(JS '4 vacuum 
with A = 0, given by 

Guv = 2r^ = , xuv = const, Tuv = const, (25) 

Auv 

while in the IR region far from the boundary, the scalar approaches the extremal value 
Air = log(2 -|- a/5) and the equations are solved by 

14 

Gm = xm = const, and Tir = 0. (26) 
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By rescaling the time coordinate one sees only xuv ~ XiR is physical, bnt for convenience in 
hnding solntions we will allow both to be free. Fixing 'hjR = 0 allows us to identify 'huv as 
proportional to the chemical potential /i of the U{l)b conserved current. 

A useful invariant of the domain wall solution is the index of refraction n, dehned by the 
ratio of the speed of light in the UV and IR CFTs: 


_ '^uv _ IR |(xiR-xuv) _ .pl(xiR-xuv) 

Rir -^uv V 7 


(27) 


This quantity is invariant under coordinate transformations, and characterizes the causal 
properties of the emergent IR conformal dynamics. 

To construct the domain wall geometries, it is convenient to consider IR-irrelevant pertur¬ 


bations about the hxed point solution (26). These perturbations can be used to numerically 
integrate away from the IR critical point at r = 0 along the radial direction, tracking the 
RG flow “upstream” to the UV hxed point at r = cx). To identify these perturbations, one 
performs a linearized huctuation analysis by substituting into the equations of motion the 
following ansatz for the IR form of the bulk helds: 


G{r) 

14/ \ 

= —r^(l + SGGj , 

(28) 

X{r) 

= XiR + Sxr^ , 

(29) 

4/(r) 

= 5T , 

(30) 

A(r) 

= log ^2 -I- a/sI -|- 5A r" . 

(31) 


The requirement that the perturbations represent irrelevant deformations to the IR hxed 
point constrains the various exponents. Specihcally, one requires a, (3,^ >0 and 7 > — 2 . 


Substituting the huctuations (28 31) into the equations of motion (21 24) one hnds that 
the linearized equations decouple, can be solved non-trivially by 


SG = 6x = 0, a=\[^-- and = (32) 

^ ’ V 28 2 ^ V 28 2 ’ ^ ^ 

and that the amplitude 54/ can be rescaled to any value under symmetries of the equations 
of motion. Thus the IR deformations are described by a single free parameter, 5A, which we 
will tune to produce domain wall solutions with various UV asymptotics. 
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The asymptotic mass of the scalar held A near the boundary is 


Id^V 


m. = - 


2 aA2 


A=0 




(33) 


which implies that near the UV boundary the scalar behaves like 

, / \ Ai Ao 

A(r -)■ cx)) ~ — + — + . 


(34) 


We will consider two particular cases, solutions in which either Ai or A 2 vanish. For reasons 
we will explain in the next subsection, we will refer to these as “Massive Boson” and “Mas¬ 
sive Fermion” backgrounds, respectively. In practice, it is straightforward to produce such 
solutions by integrating the equations of motion from very near the IR hxed point to the UV 
boundary. This involves employing the scaling symmetries of the equations of motion to hx 


and xirj then using the huctuations dehned in (28 31) to produce IR boundary condi¬ 


tions for the numerical integration of (21 24) for many choices of 6X. After each successful 


integration throughout the bulk, one can £t the near boundary behavior of the numerical 
solution obtained for A to the form given in (34), and subsequently extract the values of Ai 


and A 2 , Xuv and Tuv characterizing that solution. 

“Massive Boson” and “Massive Fermion” solutions constructed from this procedure are 
shown in hgures and ^ The Massive Boson solution, with nonzero A 2 , has interesting 
similarities to the extremal AdS Reissner-Nordstrom solution, and in some sense is “almost” 
AdSRN. Extremal AdSRN is characterized by an AdS 2 x near horizon geometry, which 
manifests as a double pole in the metric function g^r- From hgure a similar feature can 
be seen around r k. A where G very nearly vanishes quadratically, before reverting to a 
nonzero value. Moreover, the hgure shows that nearly all the scalar hair is bunched behind 
this “almost” horizon. Perhaps not surprisingly, similar properties have been observed in 
the extremal limits of various holographic superconductors studied in the literature [TS] . 

The index of refraction in this solution is large, n = 26.900, which implies that the 
effective speed of light is very slow in the IR CFT compared to the UV theory, in turn 
suggesting that the IR dynamics is nearly 0-1-1 dimensional, reminiscent of the Semi-Local 
Quantum Liquid [32] • This is another sense in which the solution is “almost” extremal 
AdSRN, since the black hole horizon corresponds to an n —)■ cx) limit. We will learn in the 
next subsection that A 2 is proportional to a dimension-2 source for a scalar bilinear; the 


^It is possible that these solutions are not unique, even up to rescalings; there may be additional solutions 
with nodes in A. If such additional solutions exist, they are probably unstable toward bosonic perturbations. 
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Figure 1: The Massive Boson background, with 5T = 1 and xiR = 4. The dashed lines in 
the plot of Gjr^ are at 14/3 and 2, indicating the values obtained in the IR and UV AdS 4 
hxed points respectively. The ratio of the speed of light in the UV CFT compared to that 
of the IR theory is n = 26.900. 
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Figure 2: The Massive Fermion background, with 54/ = 1 and xir = 4. The dashed lines in 
the plot of Gjr^ are at 14/3 and 2, indicating the values obtained in the IR and UV AdS 4 
hxed points respectively. This geometry is characterized by uuv/'f^iR = 1.861. 
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dimensionless ratio of the source to the U{1) chemical potential can be measured to be 

, 1/2 

^^ 0.0308, (35) 

indicating that indeed this solution can be thought of as a small perturbation by A on top 
of the no-scalar background, which has extremal AdSRN as a solution. 

The Massive Fermion background has A 2 = 0, and unlike the previous case, this geometry 
is not “almost” AdSRN in any sense. The function G{r) never comes close to vanishing, so 
the solution is not “close” to having a horizon. The index of refraction is substantially closer 
to unity at n = 1.861, so the speed of light does not change that dramatically between the 
UV and the IR. In this case Ai is proportional to a dimension-1 source for a fermion bilinear, 
and we have 

1.227, (36) 

so the two massive perturbations to ABJM theory, the source and the chemical potential, 
are the same to a factor of order unity. We plot the IR light cones for the two solutions next 
to a UV light cone normalized at right angles in hgure 

These backgrounds interpolate between UV and IR hxed points that are known to be 
stable in the following sense. The ultraviolet AdS 4 is stable on account of supersymmetry; 
this corresponds simply to the unitarity of ABJM theory. It is shown in that all the 
scalar fluctuations in the non-supersymmetric IR AdS 4 geometry satisfy the Breitenlohner- 
Freedman bound uni SI]. However, we do not know of a demonstration of stability of the non- 
supersymmetric AdS 4 solution against perturbations involving non-scalars; also, stability of 
the anti-de Sitter endpoints of these domain wall solutions does not by itself demonstrate the 
stability of the whole domain wall. Nonetheless, these domain wall backgrounds are the best 
candidates available for a stable holographic dual of a hnite-density state in ABJM theory. 

3.2 Holographic Interpretation 

The domain wall backgrounds constructed in the previous subsection are horizonless solutions 
to AT = 8 gauged supergravity with a non-vanishing electric potential for the gauge held. 
Thus, we broadly expect that these bulk solutions provide a holographic description of certain 
zero temperature states of ABJM theory at hnite density. The fact that a charged scalar 
is turned on in these backgrounds implies that we are studying either a deformation of the 
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Figure 3: In units where the UV light cone is 45° (dotted black), we compare the Massive 
Boson (solid blue) and Massive Fermion (dashed red) IR light cones. 

ABJM theory by the addition of a dual scalar operator, a state of the ABJM theory with 
non-vanishing expectation values for this scalar operator, or some combination of these. 
Since we have a top-down solution, we can determine the nature of the background precisely 
using the explicit mapping between the bulk helds and various single trace operators of the 
ABJM theory. 

The truncation of the maximal gauged supergravity breaks the SO{8) gauge symmetry to 
50(3) X 50(3) X 50(2). The surviving bulk helds are all singlets under the 50(3) x 50(3), 
and carry charge only under the remaining 50(2) = f/(l). The gauge held associated with 
this 0(1) is A, the active gauge held in the domain wall backgrounds constructed above. 
The gauge symmetry present in the SUGRA theory is holographically dual to the global R- 
symmetry of the ABJM theory, and thus bulk solutions with non-zero A* = T correspond to 
ABJM theory with a chemical potential /i turned on for the conserved global U{l)b current. 
The dual U{l)b current counts monopole number, and hence takes the form 

J'^~e^"^Tr(F,A + F,A), (37) 

where F and F are the held strengths for U{N) x U{N). 

Let us now connect the scalar A in our background to dual ABJM operators. As men¬ 
tioned previously, the 50(3) x 50(3)-invariant sector of A/" = 8 gauged supergravity contains 
a hypermultiplet, corresponding to two complex scalars, which can be packaged in various 
ways. The gauged supergravity naturally gives rise to 0i, 02 which are a complex 50(2) dou¬ 
blet, the real parts being parity-even scalars and the imaginary parts being pseudoscalars. 
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We can define 0i = + iPi) and <p 2 = + *-^ 2 ), and can assemble charge and 

parity eigenstates as S = Si + 182 , P = Pi + iP 2 - I2S1 also make use of the combinations 
Cl = ^(01 - ^02) = ^ (>5^ + and C 2 = ^(0i + ih) = + iP)- 

To identify the dual ABJM operators, recall the 70 scalars of the gauged supergravity 
theory live in a 35v © 35c, each of which decomposes into S'17(4) x U{l)b representations as 

35v —y 15o © IO 2 © 10_2 , 35c —y 15o © 10_2 © IO 2 , (38) 

corresponding to the operators dual to the parity-even scalars, 

15o: IO 2 : Y^^Y^^e^\ T 0_2 : Yl^Ylf-‘^\ (39) 

as well as the operators dual to the pseudoscalars, 

15o : , IO 2 : '4’{A'4’B)e ^'", IO -2 : . (40) 

Under SU{4:) —)■ SO{3) x 50(3) the 15 does not contain a singlet, while both the 10 and 
the 10 become (3,3) © (1,1). Thus the complex scalar S and pseudoscalar P living in the 
50(3) X 50(3)-invariant truncation correspond to the U(l)b-charge 2 operators 

Os = Y^Y \^^, Op = 'il^A'il’Ae^^. (41) 

The ansatz (|^ truncates the scalar sector down to one complex scalar. The truncation can 
be described in various equivalent forms: 


Cl — 0 O 01 — *02 5 — iP -H- 5i — —P 2 , Pi — S 2 , 


(42) 


and the remaining two degrees of freedom can be identihed with A and a from (|^ as 


tanh A e*" = C 2 = V2S = iV2P. 


(43) 


We note that while the Lagrangian (13) indicates the scalar has charge in a convention 


where the gauge held is dimensionless, it is convenient for us to match the natural held 
theory convention and refer to this as charge 2. Thus our background involves a simultaneous 
turning on of sources and/or expectation values for the A = 1 operator Y^Y^e^'^ and the 
A = 2 operator , with a hxed relative phase. 
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All the scalars of the supergravity theory have the asymptotic mass = —2, lying 

in the window where both the leading terms in the near-boundary expansion (34) are normal¬ 
izable deformations of AdS 4 , and correspondingly the scalars can be quantized in one of two 
ways. Supersymmetry |10l HI] requires that the pseudoscalars in the 35c have the standard 
quantization dual to an operator with A = 2, while the scalars in the 35v must have the 
alternate quantization, and be dual to operators with A = 1. For regular quantization helds, 


the mode Ai in (34) corresponds to the source, while the subleading term A 2 corresponds to 
the expectation value; this holds for our operator. For alternate quantization helds, 

Ai is the expectation value, while A 2 is the source, which holds for . Hence we hnd 

each parameter in the solution controls both a source for one operator and an expectation 
value for the otherj^ 


Ai ~ = (ReH^), (44) 

A 2 ~ </Rey2 = — (lm'0^), (45) 

where Y^ and are shorthand for Y^Y'^e^'^ and respectively. Thus for a solution 

where Ai = 0, the background corresponds to a source for the scalar bilinear, hence the 
name “Massive Boson”, as well as an expectation value for the fermion bilinear, while for 
the “Massive Fermion” solution with A 2 = 0 we have a source for a fermion bilinear and 
an expectation value for the boson bilinear. Since either KeY'^Y^e^'^ or Ira'i/jA'il’A^^^ is 
explicitly added to the ABJM held theory Lagrangian, in either case one breaks the U{1) 
symmetry explicitly. Hence these domain wall backgrounds are not precisely holographic 
superconductors, which should involve only spontaneous breaking of the symmetry. 

We note in passing that it is possible to cast these backgrounds as true holographic su¬ 
perconductors, if we pass to an alternate quantization of some of the scalars or pseudoscalars 
and hence move away from ABJM theory to a non-supersymmetric boundary theory. In a 
quantization where all the active scalars are dual to A = 2 operators, Ai = 0 backgrounds 
involve only expectation values of the dual operators; conversely A 2 = 0 backgrounds have 
no sources if all the active scalars are dual to A = 1 operators. In these two scenarios, the 
solutions are in fact holographic superconductors in the usual sense. The former case can be 
obtained as the infrared limit of a deformation of ABJM theory by a relevant double trace 
operator, essentially the square of Y^Y^e'^'^. We prefer, however, to work in the quantiza¬ 
tion dual to the ABJM theory when analyzing fermionic Green’s functions, because there 

^Additional finite boundary counterterms could in principle shift these relations at the nonlinear level, 
but the terms required by supersymmetry found in [42j vanish in our background. 
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can be no donbt abont the operators dnal to the snpergravity fermions. In this approach, 
we cannot claim to be analyzing fermionic response in a trne holographic snpercondnctor, 
bnt in a member of a broader class of symmetry-breaking backgronnds. 


3.3 Conductivities 


Before we tnrn to the task of probing the fermionic properties of these SUGRA backgronnds, 
it is sensible to wonder what lessons we can learn from the linear response of bosonic probes. 
An obvions candidate is the condnctivity of ABJM matter charged nnder the global f/(l)b. 
The imaginary part of this condnctivity appeared previonsly in [26] , where it was claimed that 
the l/cj pole in the imaginary part of the DC condnctivity was indicative of snpercondnctivity 
in the bonndary gange theory. We will briefly revisit this claim before tnrning to the real 
part of this U{1) condnctivity. 

The linear response of the gange theory current to an applied electric held is encoded in 
the retarded Green’s function, which in turn dictates the AC conductivity, a{u)): 


aiui = 


(y) 

Err. 


UJ 


(46) 


Roughly speaking, (46) shows that the real part of the AC conductivity gives a measure of 


the density of states for charged matter at zero spatial wavenumber. 

From a bulk perspective, computing this conductivity is by now a standard exercise 
in applied holography. The computation begins by turning on the coupled perturbations 
A ^ A + 6A and g ^ g + 6g where 


6 A = 6Arr{r)e 6g = 6gtx{r)e dtdx , 


(47) 


and continues by solving the equations of motion linearized about these perturbations. The 
near boundary behavior of dA^ fully determines the current-current correlator If 


A/t(i) 

6Ax{r —)■ cx)) ~ dA^^^ -\ - - —(-... 


then 


Gjj =2 

’Jx <Jx 


dA 


( 1 ) 


dA 


(0) ’ 


(48) 


and using (46) results in the conductivity. 


In hgurej^the real and imaginary parts of the AC conductivity are shown for charge trans¬ 
port in both Massive Boson and Massive Fermion backgrounds. From the rightmost plot, it 
is immediately clear that Imcr ~ l/a; at low energies. By the Kramers-Kronig relations, this 
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Figure 4: The real (left) and imaginary (right) AC conductivity in Massive Boson (darker) 
and Massive Fermion (lighter) backgrounds. The imaginary part of the conductivity has been 
multiplied by u to highlight the 1 /u pole at low energies giving rise to the delta function in 
Re a. 

necessarily implies a delta function contribution to the real part of the conductivity. Such 
a delta function does not imply that the backgrounds we are studying are holographic su¬ 
perconductors, since any translationally invariant background with non-zero charge density 
will show a similar delta function peak in the conductivity [H]. Indeed, as we saw in the 
previous subsection, our backgrounds are not true superconductors for ABJM theory, since 
the U{1) is explicitly broken, although the same conductivity calculations apply to the true 
holographic superconductors associated to non-supersymmetric alternate quantizations. 

Re (^{u) is characterized in both backgrounds by the aforementioned infinite DC contri¬ 
bution separated from the conformal plateau at high energies by a soft gap. The fact that 
the Massive Fermion background gives rise to a broader gap can be used to argue that this 
state has an enhanced suppression of charge carrying states at intermediate energies rela¬ 
tive to the Massive Boson background’s dual. As indicated above, this suppression can be 
inferred only for the states near the origin in momentum space, and thus is of limited utility 
for uncovering what is happening to the fermionic degrees of freedom in the dual ABJM 
state. This is because the natural expectation for a system of fermions at finite density is to 
organize into a Fermi surface at some finite k = kp. Thus, to address questions related to 
the fermionic nature of these states it would be more appropriate to study current-current 
correlators at non-zero k along the lines of [43l 03]. An alternative approach, which we 
adopt in the present work, is to study the fermion response of the ABJM states directly, 
using appropriate fermion probes. 
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4 Fermion Response in the Domain Wall Solutions 


We would like to study the linear response of states in the gange theory to the insertion 
of varions fermionic operators, which is characterized by an assortment of fermionic two- 
point fnnctions. Holographically, these two-point fnnctions are compnted from the linearized 
flnctnations of snpergravity fermions abont the classical (bosonic) backgronnds of interest. 
Here we will describe the calculation of these eqnations, the mixing between modes, and will 
discuss a simplihed fermion mixing matrix. 

4.1 Coupled Dirac equations and holographic operator map 

We will focns on spin-1/2 helds that cannot mix with the gravitino sector. Under the 
SO(8) —)■ SU(4) X I7(l)b SO(3) x S'0(3) x SO(2) decomposition, we have for the gravitini 
in the 8s, 


8s —t 6o © I2 © 1-2 (3,1)0 © (1, 3)0 © (1,1)2 © (1, l)-2 , (49) 

and thus we can avoid mixing in the SO(3) x S'0(3)-invariant backgronnds as long as we 
stndy fermions in representations other than (3,1), (1,3) or (1,1). The spin-1/2 helds live 
in the 56s, which decomposes as 

56s —t 152 © 15-2 © IOq © IOq © 60 

(8, 3)2 © (3,1)2 © (1, 3)2 © (3, 3)_2 © (3, l)-2 © (1; 3)_2 © (50) 

2(3,3)o©2(l,l)o©(3,l)o©(l,3)o, 

and thus we see there are fermions in the (3,3) of SO(3) x SO{3) that cannot mix with 
the gravitini. They may, however, mix with each other, and in fact do, as we shall see. 
The different SO{2) charges of the fermions in the (3,3) representations are no obstacle to 
mixing because the SO{2) symmetry is broken by non-zero A in onr backgronnds. These 
(3,3) fermions are not in the fermionic sector of the SO{3) x SO{3) trnncation discussed 
previously, and thns to obtain their dynamics we must return to the fnll A/” = 8 snpergravity 
theory. 

Dropping gravitino pieces that will not conple, the relevant terms for the spinor Xijk = 
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X[ijk] in the full M = 8 gauged SUGRA Lagrangian are |HU] : 

^+ hx.) 

+ + h.c.) ^ (51) 

The fermion tensor is defined through 

. (52) 

and the covariant derivative is 


DfxXijk — '^iiXijk „ i Xljk 


Ib.‘ 


f_i j "X^ilk 




^ k 


(53) 


where is the covariant derivative dehned with respect to the spin connection u^at'- 


v^ = d^--u,^r“‘, 


(54) 


and the composite SU (8) connection B is determined through the vanishing of the diagonal 
blocks of Q to be 

13^) = ^ - v^^^^%VjkLM^ , (55) 

where is covariant only with respect to SO{8), and does not act on SU{8) indices i, j. 

We now consider the Dirac equations for the Xijk in the (3, 3). If we take the hrst 50(3) 
to act on greek indices a = 3,4,5, the second 50(3) to act on roman indices a = 6,7,8, 
and the 50(2) which corresponds to the active gauge held to act on hatted indices a = 1, 2, 


the fermions that transform as four distinct copies of (3,3) in (50) are readily seen to be 


those of the form Xapc Xabc, and Xai 3 c, where we recall the 50(3) antisymmetric product 
3 3 = 3. Thus, an example of a set of fermions with the same 50(3) x 50(3) quantum 

numbers is {x 467 , Xsss; X 4 i 8 ) X 428 }; these helds may all mix with each other, but not with any 
others. We will study these four fermions; any other analogous quartet has results related 
by group theory. We assemble the fermions into charge eigenstates, 

X2 = X428 + W418, X2 = X428 — W418, Xo = X467 + W538, Xo = X467 — fX538 • (56) 


These fermion helds (56) are dual to spinor ABJM operators of the form with A = 3/2 


and the 56s arising in the product 8v x 8c = 56s © Sg. Under 5f/(4) x f/(l)b the 56s 
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decomposes as 


56s —^ 152 © 15—2 © lOo © IOq © 6o , 


(57) 


corresponding to the ABJM operators 


152 : 




. 2 t 


15_ 


2 • 



(58) 


lOo: , lOo: , 60 : comb, of and Y^^ipsi , (59) 

where the other linear combination of the antisymmetric part is part of the Sg. Under 
SU{A) SO{3) X ^0(3), the 15, 10 and 10 all contain a (3,3). The helds X 2 and X 2 are 
charged under 17(1);,, and hence will he in the 152 and 15_2, while Xo and Xo sit in the IOq 
and IOq. Tracing through the indices one hnds 


X 2 


(UV2 - UVi + - UVa) e"" , 

(60) 

X 2 



(61) 

Xo 

y-)- 

, 

(62) 

Xo 


U1V4 + nVi - U2V3 - U3V2 . 

(63) 


The Dirac equation for these fermions can be obtained from (51), plugging in the values for 
the supergravity quantities described in section appropriate to the backgrounds discussed 
in section IH The result takes the form 


(*r'^v^i + s)x = o, (64) 

where 1 is the identity, y is a 4-component vector containing the spinors, and S = A-l-P + M 
with A, P, and M describing gauge, Pauli, and mass type couplings, respectively. We hnd 
that these matrices fail to commute, and the four spinors mix nontrivially. In this paper, we 
study a simplihed version of the fermion interactions, which we call the non-chiral mixing 
matrix, as a step towards the full interactions. These interactions capture many aspects of 
the full top-down mixing, but differ most signihcantly in containing no chiral Ts terms. The 
full, chiral interactions will be considered in [28] 

^An earlier version of this paper mistakenly identified this as the full mixing matrix. 
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The non-chiral matrix S is 


^ —1^(3-|-cosh 2A) 

0 — sinh A 

0 

\ 

0 

^^(3-|-cosh2A) 0 

— sinh A 


— sinh A 


^^sinh^ A 


V 0 

— sinh A sinh^ A 

i jt 

/ 


where we have written ^ Jf = One can see that X 2 and X 2 are charged, 

and as with the charged scalar, we identify this as charge ±2 in the natural normalization 
of the held theory. Meanwhile xo and xo are neutral, but have Pauli couplings to the held 
strength. This basis is diagonal at the ultraviolet hxed point A = 0, where all the fermions 
are massless. As the scalar turns on away from the boundary, it rescales the gauge couplings 
of X 2 and X 2 and, most importantly, introduces couplings between diherent fermions: there 
is an interaction between the neutral and the charged fermion of the schematic form (j)X 2 Xo 
with 0 the charged scalar, and a coupling between the neutral fermion and its conjugate of 
the form A|0py;oXo- 

It is interesting to compare our Dirac system to other fermionic equations used in holo¬ 
graphic superconductors. In [12] an ordinary Dirac equation with tunable charge and mass 
was studied, and the superconducting character was inherited from interactions with the 
background. In [20], new terms were added to the Dirac equation to emulate the effects of 
the Cooper pair condensate by coupling the spinor to its conjugate, with “Majorana” terms 

AC^ (j)*X^C{r] + r]5T5)x + 'i^-c., ( 66 ) 

leading to a Dirac equation of the form 


(ir'V^ - m + gT^A^) x + (>) + rkT^)<l)Bx' = 0, 


(67) 


where rj and rj^ are coupling constants, Ts is the chirality matrix, and B is related to the 
charge conjugation matrix via C = The scalar must have = 2g^, and its con¬ 

densation breaks the U{1). This leads to terms like XX + X*X* fh® effective Lagrangian, 
analogous to the cc -|- A A terms in a BCS superconductor Lagrangian. 


Our system can be viewed as an elaboration of (67). Instead of coupling a single charged 


held to its conjugate, our system has a “Cooper pair” coupling between the charged scalar, 
the charged held X 2 and the neutral held xo, breaking gauge invariance when the scalar con¬ 
denses, as well as a “Majorana” coupling between the neutral held xo and its own conjugate 
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Xo, mediated by the gauge field and the scalar squared. Looking at a coupling like (66), it 


is somewhat natural to think of the scalar field as being dual to the Cooper pair fermion 
bilinear. Our system is a little more complicated: the two fermionic operators are of the 
form Y'lp, so the Cooper pair is some part of Y'lpY'ip, while the operator dual to the scalar 


that condenses is either of the form Y^ or -0^. The non-chiral mixing matrix (4.5) lacks the 
Ts terms found in [201, but these are present in the full chiral mixing matrix. 


In the next subsection we discuss solving the Dirac equations (64, 4.5) and relating the 


results to Green’s functions for these operators, from which we can calculate the normal 
mode spectrum and the spectral functions. 


4.2 Solving the Dirac equations and spinor Green’s functions 

The analysis of Dirac equations in nonzero density backgrounds is by now standard in the 
literature, for more details see for example IHEIEIE]. We rescale the spinors by a factoij^ 
to cancel the spin connection term in the Dirac equations, and Fourier trans¬ 
form as with frequency oj and spatial momentum k chosen to lie in the x-direction. 

Next, we make a convenient choice of Clifford basis where the relevant T-matrices are block 
diagonal, 

^(Tl 0 
0 (Ji 


pr _ 



r* = 


px _ 



( 68 ) 


n„ = - 
2 


By then defining the projectors 

1 - 

we can write the four components of the spinor y as 

Xoit — nQ-Tj-y, 


p^^-(i±ir) , 


(69) 


(70) 


with a = 1,2. With our choice of Clifford basis, it is fairly easy to see that the Dirac 
equations do not mix spinor components with different a, meaning we can split them up into 
two decoupled sets of equations, and one can show the solutions of the two sets are related 
simply by fc —)■ —k. We also note that our Dirac equations have a discrete conjugation 
symmetry, being unchanged under the simultaneous substitutions 


uj —y —UJ , 


k —)• —k , 


X+ ^ X+ , 


x_ O -x_ , 


(71) 


^The metric function x appearing in this factor should not be confused with the spinor. 
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where x± represents the ±-components of both X 2 and Xo- Thus we can restrict to A; > 0 
and a = 1 (dropping the a-label), using 0 to reconstruct A; < 0 and obtaining a = 2 
simply by changing the sign of k. 

Even having restricted to half the spinor components, our system still involves eight cou¬ 
pled hrst-order equations. Beginning in the deep IR, as per the usual holographic dictionary 
we want to impose appropriate boundary conditions to compute retarded Green’s functions. 
In the infrared limit, the coupling matrix S is off-diagonal in the charge basis {x 2 , X 2 , Xo, Xo}, 
but becomes diagonal in the “mass basis” 


XW — X538 — X416, XX — X467 ~ X426, XY — X538 + X416, Xz — X467 + X426 , (72) 


which diagonalizes the mass matrix M. As r —)• 0 each of the mass basis spinors a = 
W, X, Y, Z obeys a second order uncoupled equation of motion of the form 


0 


Xa + ~X'a + 
r 


p. 


Lm ’^ir(1 

A ' 


^ir) \ 




Xa 


(73) 


where miR is the dimensionless fermion mass at the IR hxed point: 


miR = mLiR = ±\l 


(74) 


and p is the 4-momentum combining u and A: in a way respecting the IR Lorentz invariance: 




u 


(75) 


-'IR 


The character of the solutions depends strongly on whether p is timelike or spacelike. When 


p is spacelike, (73) admits solutions that are either regular or divergent in the IR. The regular 


solutions are of the form 


Xa{r) = Na—^K^l_ 
A/r 


miR 




(76) 


with Na a normalization constant. On the other hand, for timelike p, the solutions are 
oscillatory in the IR, being either infalling or outgoing. 

In the far UV, the charge basis spinors decouple and solve massless second order equations 
of the form 


Xct + ~Xa — 0 , 
r 


(77) 
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where the index a now stands for the i-components of each charge basis spinor, and we 
are suppressing an index labeling the distinct elements of the charge basis. The leading 
constant solution^ for the X- modes are associated to the expectation values {O) of the 
dual operators, and those for the x+ modes are associated to the sources J: 

X+(r) ~ Jiu, k) + 0{r-^) , x-{r) ~ (0(a;, k)) + 0{r-^). (78) 


The choice of which of x± is associated with the source and which with the expectation value 
is determined for ABJM theory by supersymmetry uni sn Dro]. 

Were the fermions decoupled, we could solve the Dirac equation for just one of them with 
the others vanishing; imposing suitable boundary conditions in the IR would compute the 
relationship between that dual operator’s source and its expectation value. In our system 
this is not the case; a general solution to the system of equations leads to all four sources 
J and all four expectation values {O) turning on. Considering the response of the four 
expectation values to varying the four sources, we we obtain a matrix of Green’s functions, 
which schematically takes the form 



5{p^) 

6J^ 


j*=o 


(79) 


To properly define this matrix Green’s function, we follow a recipe very similar to the one 
advocated in [25l HS], searching for solutions to the equations of motion with suitable IR 
boundary conditions in which only one ABJM operator is sourced at a time. Given such a 
solution, standard application of the AdS/GFT dictionary for spinors allows us to read off 
the linear response of the operators to this source, and the associated entries in the matrix 
Green’s function. 


In practice we proceed as follows. The IR normalization constants Na (76) can be chosen 
independently for each of the bulk spinors. This is guaranteed by the linearity of the equa¬ 
tions of motion combined with the fact that the bulk fermions completely decouple in the 
IR. Imposing the proper boundary conditions in the IR, one can vary the Na and see how 
the sources J in the UV change. In this way we can construct a linear map T between the 
IR data N and the UV sources J = (J"^, J^, J^): 


TN = J. 


®The rescaling described above (68) removed a leading factor of r 


(80) 
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Figure 5: Cartoon of regions with dispersionless modes, meaning regular, rather than in¬ 
falling, boundary conditions in the IR. At left is an extremal horizon, such as AdSRN; in 
the middle, an IR AdS region, as in the present paper; at right, an IR singularity. 


The inverse of this map allows us to construct the IR data needed to produce a bulk solution 
with any desired values for the dual sources. Once such a solution is known, sources and 


expectation values can be read off using (78) and plugged into (79) to obtain the Green’s 
function matrix. From a practical standpoint, constructing the 4x4 matrix T is a straight¬ 
forward but computationally tedious affair. One can completely determine the 16 complex 
entries by integrating the equations of motion four times, with four distinct (but arbitrary) 
N. After each integration the values of J are computed from the UV asymptotics of the 
solution, eventually yielding 16 equations for the unknown entries of T. This process must 
be repeated for each value of {u, k) of interest. 

The appropriate boundary condition in the IR depends on whether the IR 4-momentum 


p (75) is timelike or spacelike. If timelike, the choice of infalling boundary conditions leads 
to calculating retarded Green’s functions. This boundary condition is complex, leading to 
a non-Hermitian matrix of Green’s functions. Solutions that vanish at the boundary are 
quasinormal modes and are associated with poles in the Green’s functions at complex u and 
corresponding excitations with hnite lifetime. This occurs “inside” the lightcone in the u-k 
plane. 

If p is spacelike, the “infalling” boundary condition used to compute retarded correlators 


can be analytically continued to the regular solution (76). This is a real boundary condition 
on the mode in the IR, which leads to a Hermitian matrix of Green’s functions. Solutions 
that vanish in the UV as well are normal modes, and are associated with poles in the Green’s 
functions at real u, and excitations that are dispersionless. 


Thus, the light-cone structure given by (75) divides the u-k plane into a region inside the 
light-cone where modes decay, and a “stable wedge” outside the light-cone with dispersionless 
excitations. We will see this in the next subsection, where we will focus on the normal 
modes and associated dispersionless excitations. This can be contrasted with other types of 
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geometries: Reissner-Nordstrom and its cousins with regular horizons have the light cone hll 
up the entire k-u plane, and hence have unstable modes except potentially at a; = 0 itself, 
while the IR singular geometries of |9], HH] have a dispersionless region for |a;| < A for a 
constant A, independent of fc; these are contrasted in £gure|^ 


4.3 Fermion Normal Modes 

The matrix T{u},k) contains all of the information necessary for identifying the locations 
of any fermion normal modes which may appear in the bulk. Any such normal mode can 
be dehned as a solution to the equations of motion which decays in the far IR and whose 
“source” falloff in the UV vanishes—in other words, a regular solution to the bulk spinor 
equations at some (cjat, k^) such that J = 0. Clearly any non-trivial solution with this 
property implies a zero eigenvalue of T, and thus one discovers that 

det T(a;Ar, ^at) = 0 . (81) 


This expression provides a powerful method for locating the fermion normal modes, and can 
be used to determine their location to very high accuracy. 


The result of applying the diagnostic tool (81) to the fermionic perturbations of the 


Massive Boson domain wall solution with the non-chiral mixing matrix (4.5) is shown in 
£gure[^ and for the Massive Fermion solution in figure We plot the results for k > 0, but 
due to (71) the spectrum is invariant under {k,u) —>■ {—k, —u). As anticipated, the normal 
modes appear in bands that are conhned to the exterior of the IR lightcone, inside the “stable 
wedge”. Our numerical search reveals two bands for each domain wall solution within the 
kinematic regions showiil^ In both cases, one of the bands passes through a; = 0, indicating 
that the non-chirally coupled fermions include gapless fermionic degrees of freedom. Because 
these gapless fermionic modes appear at hnite momentum, their presence indicates that the 
dual fermions organize into a Fermi surface. These ungapped bands in both cases begin at 
the upper boundary of the light cone, and asymptote along the lower edge of the light cone 
as far as our numerics can follow. Both cases also possess a gapped band, which appears 
to both begin and end along the lower light cone edge. In the Massive Fermion case, the 
gapped and ungapped bands come close to each other along this edge, but the gapped band 
appears to terminate before they coincide. 


®Wliile we do not completely exclude bands of normal modes at higher momentum than what is shown 
in the figures, a rough numerical search along the edges of the IR light cone revealed no further interesting 
features for kvuv/^J- < 10- 
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Figure 6: The band structure of fermion normal modes in the Massive Boson background. 
The shaded blue triangle is the stable wedge where it is possible for normal modes to appear, 
and the solid blue curves are the locations of fermion normal modes of the bulk theory, as 
determined by solving (81). The intersection of the dashed line with one band indicates 
the presence of a gapless mode. This band appears to terminate where it reaches the top 
boundary of the shaded region, but follows it closely along the bottom edge as far as our 
numerics allow us to compute. 


Thus we see suggestions of both gapless and gapped excitations of ABJM collective 
fermionic degrees of freedom, which since they correspond to poles at real u are perfectly 
non-dissipative, at least at large N. We have set the scales of both hgures so that details can 
be seen, but it should be remembered that the wedge outside the lightcone where such stable 
fermionic excitations can exist is much smaller for the Massive Boson case than the Massive 
Fermion case (compare the light cones in hgurej^. Inside the IR lightcone no additional 
normal modes exist, but only quasinormal modes at complex frequencies corresponding to 
excitations that decay; we will get a sense of such modes in the next subsection. 

As discussed in section there is a correlation between the strength of the symmetry 
breaking source and the size of the IR light cone. When the U{1) symmetry breaking is 
turned off entirely and only the chemical potential remains, the solution is AdSRN and the 
“IR light cone” effectively hlls the oj-k plane, leaving no space for stable modes; indeed, as 
hrst demonstrated in |1] and shown for top-down ABJM fermionic fluctuations in [7], this 
geometry supports fermionic zero-energy modes at hnite momentum, but no stable (inhnitely 
long-lived) excitations at hnite frequency. As the symmetry breaking is turned on weakly in 
the Massive Boson case, the light cone closes slightly and a kinematic wedge appears where 
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Figure 7: The band structure of fermion normal modes in the Massive Fermion background, 
determined by solving (81). Again, there is a gapless mode at hnite momentum. At higher 
momentum the gapped band approaches the ungapped band, but appears to meet the IR 
lightcone before the two bands coincide. As in the Massive Boson background, the ungapped 
band traces the bottom edge of the stable region as far as our numerics can reliably follow 
it. 


stable modes exist; the Massive Fermion case has symmetry breaking of the same order as 
the chemical potential and a much larger stable wedge. 

It is tempting to conclude that in turning on the symmetry breaking source, some sector 
of the gauge theory mediating decays of the fermionic excitations has become gapped. That 
the gap is dehned by the boundaries of the IR lightcone and not by the size of the symmetry 
breaking deformation alone can be understood as a consequence of the emergent IR conformal 
symmetry. Far in the IR, the only relevant dimensionless scale is dehned by the huctuation, 
like Air = u/kvi^. Dialing up the strength of the symmetry breaking source closes the 
IR lightcone further. While the symmetry breaking source is not the same operator in 
the Massive Boson and Massive Fermion geometries, being a scalar bilinear with monopole 
operators in one case and a fermionic bilinear with monopole operators in the other, we may 
speculate that in these geometries it is the size of the symmetry breaking rather than the 
details of its nature that most strongly inhuences the IR dynamics. This can be inferred 
from the fact that both sorts of deformations drive the UV theory to the same IR hxed point. 

Geometries with “good” IR singularities were studied in [9] , and in those cases inhnitely 
long lived fermionic excitations also appeared, again sometimes connected to a zero-energy 
mode at hnite momentum. In that case the stable region was not a wedge, but a band 
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defined by |a;| < A. Beyond the value A (which is proportional to the chemical potential) 
the normal modes were found to move off the real u axis and the fluctuations consequently 
acquired a hnite width. It is likely something similar happens in the present case as well. 

The normal mode analysis does not, in and of itself, provide any information about 
which ABJM fermions are participating in these excitations. By virtue of our top-down 
holographic approach to this system, we can address this question, and at the same time 
better understand the fate of the normal modes beyond the boundary of the stable region. 


4.4 Spectral Functions 


The calculation of the normal modes (81) treats all four fermions symmetrically. However, 
the fermions do not all participate in each mode equally. A normal mode may be thought 
of as a solution for which all the sources vanish; however, the four expectation values may 
behave differently, as some may vanish in the normal mode and some may not. Equivalently, 
one may imagine approaching a normal mode in the u-k plane while keeping a source hxed, 
and some expectation values will then diverge. Expectation values that are nonzero in a 
normal mode will thus be associated to poles in the matrix of Green’s functions. It is 
interesting to determine which fermionic operators participate in which collective normal 
modes. 

A natural way to explore this is to study the spectral function matrix, proportional 
to the anti-Hermitian part of the retarded matrix Green’s function i{G{oj,k) — G^{u,k)). 
This spectral function quantihes the fermionic degrees of freedom at a given frequency and 
momentum which overlap with the fermionic operators of the ABJM theory. Unlike the 
normal mode bands found in the last subsection, the spectral function will be nonzero outside 
the stable wedge, and will provide a sense of the existence of unstable modes in this region. 
However, inside the stable wedge the spectral function itself is hard to examine. This is 
because since all excitations there are perfectly stable, the spectral function is zero except 
for delta function singularities; our numerical solutions cannot pick up these peaks, meaning 
the plots of these regions would be quite boring. To remedy this problem, we recall the 
Kramers-Kronig relations require the real parts of the Green’s function matrix to possess 
l/cj-type poles when the imaginary parts have delta functions. Hence we choose to study 
the quantity G'^G{u}, k), which will bring together both the real and imaginary parts of the 
Green’s functions. We can then plot trG^G, which will be a basis independent quantity 
capturing both excitations in the stable wedge outside the IR lightcone and hnite-width 
excitations inside the IR lightcone. 
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Finally, we can define matrices that project onto the subspaces of dehnite charge: 


P+= diag{l, 0,0,0} P_ = diag{0,1,0,0} Pq = diag{0, 0,1,1} 


(82) 


(c.f. equation (56)). Then, irP+G^G will measure the excitations of X 2 alone, etc. 

In hgures and various projections of tr for the non-chirally mixed fermions are 
plotted for the Massive Boson and Massive Fermion backgrounds, respectively. Our hrst 
observation is that the spectral density inside the stable region is strong along the curves 
of the normal modes, as we would expect. This density continues outside the stable wedge, 
pointing to the presence of nearby unstable modes. Most of the bands of density along 
the normal mode curves are strong; the exception is the region of the gapped band in the 
Massive Fermion background, which is weak enough to show up in our plot as a series of 
distinct points. We continue to plot A: > 0, but the spectrum is again invariant under 
{u},k) —>■ (—cj, —k)] thus the band in the lower-left corner of the Massive Fermion plots is 
the continuation of the band that exits the plot in the upper left. 

By projecting onto the charged and neutral subspaces, we can identify which fermions 
participate in which bands of spectral density and associated normal modes. For both 
backgrounds, it is X 2 and xo fhat dominate both the gapped and gapless bands of normal 
modes for A: > 0. X 2 (and xo) participates only very slightly along these same curves. We 
stress that this asymmetry is in part arbitrary: the conjugation symmetry (0 tells us 
that X 2 and Xa will have similar strong excitations for (w, k) —>■ (—w, — A:), or in the a = 2 
components with u —)■ —uj. We note that for the Massive Boson background the gapped 
band continues in the P_ projection out of the stable wedge and into the light cone. The 
symmetry described above also involves charged conjugation in the charged subspaces, and 
this curve can be seen to hnish in a small tail just above the vertex of the light cone in 
the P+ projection; this feature will be shown to be a remnant of a stronger band when we 
consider modifying the couplings in the next subsection. 

It is interesting to quantify the relative participation of different fermionic modes at a 
particular point on these bands; we choose to look at the Fermi surface point along the corre¬ 
sponding curve, at zero frequency (relative to the chemical potential) but hnite momentum 
k = kp- One can turn on a unit source at cj = 0 and k = kp ior each of the four supergravity 
fermions and catalog the response of the system to this source in the rows and columns of 
the retarded Green’s function. Diagonalizing the Green’s function at a; = 0 and k = kp 
explicitly reveals an eigenmode with diverging eigenvalue. Denoting this eigenmode ^kp, one 
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Figure 8: Plots of G^G for the Massive Boson background. Within the wedge marked by red 
edges, all excitations are stable. 
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Figure 9: Plots of for the Massive Fermion background. Within the wedge marked by 
red edges, all excitations are stable. 
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Figure 10: The squared amplitude |c/p of each ABJM operator participating in the “Fermi 
surface” zero-energy mode in the AdS-RN background with no scalar (left), and the Massive 
Boson (middle) and Massive Fermion (right) backgrounds. The normalized Fermi momentum 
kpvuy/fi in the three cases are 0.53 (AdS-RN), 0.58 (MB), and 0.48 (MF). Note that in all 
cases the contributions from are insignihcant. 


can then write 

^kp = '^cix^ (83) 

I 

where the the bulk supergravity fermions dual to the ABJM operators we study. The 

amplitudes c/ thus quantify the amount in which various ABJM fermions are involved in the 
fermionic zero-energy mode. This decomposition is shown in hgure [TO] for the Massive Boson 
and Massive Fermion backgrounds, quantifying how the normal mode is primarily composed 
of X 2 and Xo] this is readily understandable as being a result of the direct mixing of these 
modes due to the symmetry-breaking “Cooper pair” coupling of the form 0XoX2- We see 
also that xo and X 2 barely participate at all. This suggests that the Majorana self-coupling 
term |0PxoXo does not have much effect on the collective mode at the Fermi surface. 


4.5 Modifying couplings 

It is interesting to ask how these results change as we modify the background or the couplings. 
This can give us an idea of which couplings are “responsible” for the effects we see. For 


example, hgure suggests that the charged-neutral coupling is much more important than 
the neutral-neutral coupling, because the X 2 -X 0 niixing is strong and the Xo-Xo niixing is 
weak, and we will see that indeed this is true. It should be kept in mind that all these 
modihed couplings, including the non-chiral mixing matrix we have primarily studied, take 
us outside the top-down approach, since we do not know any explicit embedding in M- 
theory for these fermion equations. We will comment on three modihcations of our system, 
as follows: 

1. We will keep the Lagrangian the same (i.e. still use W = 8 supergravity) but turn off 
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the scalar field A while keeping the chemical potential. The result is the AdSRN black 
brane corresponding to ABJM theory at zero temperature, deformed only by chemical 
potentials. In the terminology of [10], this is the four-charge black holej^ 

2. We will consider massless charged Dirac fermions with no couplings to the scalar fields 
in the Massive Boson domain wall background. This is analogous to the approach of 

ra- 


3. Also in the Massive Boson domain wall background, we will modify the equations of 


motion (64)-(4.5) for A/" = 8 fermions in only one regard, namely by omitting the off- 
diagonal “Cooper pair” and “Majorana” couplings. These couplings are similar to the 
ones considered in m- 


Results for the Massive Fermion background are similar. All these interactions can be 
regarded as steps on the road from the simplest fermion couplings in the simplest finite- 
density backgrounds to the full chirally mixed interactions. 

The existence of the “stable wedge” is a property of domain wall backgrounds with IR 
AdS regions. If we pass to the AdSRN background, then the stable wedge is closed, and 
all excitations away from a; = 0 are dissipative. At a; = 0, there is a Fermi surface for the 
charged fermion [7|, and the neutral fermion is at a special transition point between a pole 
in its Green’s function and a zero as other chemical potentials of the system are varied [10] . 
We include the observation that X 2 is entirely responsible for this Fermi surface singularity 
in the AdSRN background in figure [To] 

Turning the scalar back on leads to the backgrounds studied in this paper and opens 
up a stable wedge. We expect there to be stable modes in this wedge for generic charged 
fermions. Indeed, in figure 11, we see similar lines of poles for charged, massless fermions in 
our Massive Boson background with elementary Dirac equations; only the neutral case does 
not acquire a band of stable excitations. Thus we conclude the existence of stable fermionic 
modes is a generic property of the background once the IR AdS region exists and the stable 
wedge appears. 

The Dirac equations associated to the non-chiral mixing matrix are substantially more 
complicated than these, involving additionally X 2 X 0 and XoXo couplings, as well as Pauli 
terms and a running of the gauge couplings with the scalar. Faulkner et ah [20] also dis¬ 
cussed how the turning on of a Yukawa coupling caused bands of excitations that crossed 

^Note that due to a triality rotation carried out in [10], the sum of all four gauge fields there corresponds 
to our single gauge field here. 
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Figure 11: Plots of the modulus squared of the Green’s function for massless probe fermions 
of various charges. 


to repel each other, leading to a gap in the dynamics. In that case, there was only a single 
charged fermion, and the coupling had a Majorana character coupling the fermion to its own 


conjugate as in (67). In general if the particle has a pole at momentum k, the antiparticle 
will have this pole at —k. However, one can see that the Ps factor mixes the a = 1 and 
a = 2 components, which introduces an additional flip of the sign of k; for this reason the 
authors of m preferred the P 5 interaction, which couples two modes with poles at the same 
momentum and leads to eigenvalue repulsion generating a gap. One may ask whether a 
similar principle of repulsion between bands brought on by a mutual coupling applies in our 
case. 


In hgure 12, we plot the normal mode structure for our fermions in the Massive Boson 


background, with the off-diagonal couplings in the mixing matrix (4.5) removed but the 
diagonal terms preserved, and compare it to the full top-down results previously shown 
in hgure In the left plot, describing the decoupled case, there are three bands: the 
gapless, yellow band stretching from top to bottom is associated to X 2 , while the red band 
crossing this coming from the lower edge of the stable wedge to the left is the neutral 
fermion Xo- Meanwhile there is a third band in orange, associated to the oppositely charged 
X 2 i crossing the X 2 band below the upper boundary of the wedge and the xo band close to 
the origin in u-k space. This orange band is also gapless, displaying a zero-energy mode 
around k uuv//4 = 0.09. We note the resemblance between the X 2 and X 2 bands shown there, 
and the free fermion q = 2 and q = —2 cases shown in hgure 

By comparing the band structure with and without off-diagonal couplings we can get an 
idea of how the couplings modify the bands. The lower-right crossing of bands results in 
both mixing and repulsion, as the crossed X 2 and Xo bands transform into uncrossed bands 
involving a mixing of both fermions. This repulsion, however, does not create a gap; unlike 
the simpler case in [ 20 ] there is no reason for the repulsed crossing to exist at cu = 0 , since 
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k k 

—Vuv ~ ^uv 

Figure 12: Normal mode structure in the Massive Boson background with the off-diagonal 
couplings turned off (left), compared to the full top-down result of hgure (right). With 
the off-diagonal couplings turned off there are three bands, associated to X 2 (orange), X 2 
(yellow), and xo (red), intersecting in three places. Turning on the couplings between the 
different fermions induces repulsion between the bands, as described in the text. 


it involves the coupling of two distinct fermions instead of a fermion to itself. Meanwhile, 
the two crossings involving the orange X 2 band lead to repulsion without mixing; the band 
of normal modes associated to X 2 is pushed off beyond the stable wedge, ending up as the 
small tail visible just above the light cone vertex in the P+ projection of hgure while the 
X 2 band only acquires a tiny X 2 component. As the X 2 band is pushed off in this way, its 
associated zero-energy mode disappears, thus gapping X 2 (for A: > 0). In this case, as in [20] . 
it happens that the coupling has created a gap. 

The tiny amount of mixing between the X 2 , Xo sector and the X 2 , Xo sector suggests that 
the Majorana Xo'Xo coupling is relatively unimportant, and on the whole this proves to be 
the case; turning it off alone removes the small contribution of the P+ sector to the normal 
modes, but does not change any of the overall structure. The X 2 -X 0 “Cooper pair” coupling 
is the dominant interaction. 
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5 Discussion 


Perhaps the most powerful aspect of our approach is the explicit holographic map between 
the supergravity modes in our gravitational solutions and various operators in the ABJM 
theory. This “top-down” application of gauge/gravity duality opens the door to various 
interpretations of our results in the context of zero-temperature states of a held theory 
whose operator content is well understood. The natural next step is to examine the full 
top-down chiral mixing matrix, including the Ts terms. The non-chiral mixing has already 
gapped out some of the lines of normal modes; a natural question is whether the full chiral 
mixing will gap out the fermionic huctuations completely. This question is answered in the 
affirmative in |28j . 

In our setup, we have deformed the ABJM theory in two ways which our analysis makes 
precise. The hrst is by the addition of a chemical potential for the U{l)b current, placing 
charged ABJM matter at hnite density and sourcing a relevant deformation away from the 
UV hxed point. In the standard presentation of ABJM theory, the scalars and fermions 
"ipA are neutral under U{l)b-, which is carried only by the monopole operators Accordingly, 
the composite monopole-fermion/scalar operators of ( |l7| ) hlling out the 8v,c carry U{l)b 
charge, and the natural interpretation of the states we study is as zero temperature phases 
of composite matter at hnite density. 

In fact zero-temperature phases of such composite matter have arisen in other hnite- 
density investigations of 2-1-1 dimensional held theories, beyond holography in the large-A^ 
limit. This is perhaps most famously apparent in the context of the fractional quantum-Hall 
ehect |l6], but related phases have also appeared more recently in e.g. HZlllHlllg]. Particle- 
vortex duality in three dimensions exchanges objects charged under an “ordinary” symmetry 


with those charged under a current associated to the dualized gauge held of the form (37), 
conserved by virtue of the Bianchi identity. Many theories which permit such a duality are 
more amenable to calculation in terms of the “magnetic” variables which generate J;,, and 
thus these variables can often provide a relatively simple description of complicated phases 
of strongly coupled matter. 

Our calculation of the spectral functions for composite fermions can help better under¬ 
stand the nature of these putative phases. One of our main results is the appearance of 
delta function singularities in the spectral functions within a particular kinematic window 
controlled by the properties of the IR hxed point. These hnite momentum singularities sig¬ 
nal the presence of stable excitations which overlap with the fermionic operators written 


in (60 63). One plausible explanation of these spectral features is that the hnite density of 
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composite fermions (those transforming in the 8c) have arranged themselves into a Fermi 
surface at a; = 0 and k = kp, and the IR excitations around this Fermi surface are weakly 
interacting and thus long-lived. In this picture, the low energy features of these states are 
qualitatively similar to a Fermi liquid of composite fermions]^ 

The other deformation we have dealt to the ABJM theory is the addition of a source 
which explicitly breaks the global U{l)b- In the states that we focus on, this breaking results 
in a non-vanishing expectation value for composite boson or fermion bilinears. In a sense 


developed in some detail in section 4^, it is the breaking of the monopole number density 
that permits stable excitations in the vicinity of the Fermi surface. It is interesting that 
the fermionic response indicates that the system remains gapless even though the U{l)b has 
broken. 

We are now in the position to ask how our results compare to other zero-temperature 
states of composite matter. One particularly interesting example is A/” = 4 supersymmetric 
QED, which is acted on by mirror symmetry and hence like ABJM theory permits a descrip¬ 
tion in terms of magnetic (composite) variables. In |18] it was shown that the IR physics of 
this theory with a uniform density of magnetic impurities is described by phases in which 
an “emergent Fermi surface” consisting of composite fermions organizes into a Fermi liquid. 
Moreover, it was found that this theory permits a phase in which composite bosons acquire 
an expectation value, yet the Fermi surface persists. Understanding to what extent our 
ABJM states match the expectations for these novel phases would be interesting. 

In m, it was argued that the pattern of holographic Fermi surfaces in symmetry¬ 
preserving backgrounds of W = 4 super-Yang-Mills theory and of ABJM theory can be 
predicted based on the form of the dual held theory operators. In particular, the held theory 
scalars involved in the dual operators may have expectation values, and if they do, then the 
“boson rule” of m predicts the existence of a Fermi surface. A slightly subtle point is that 
the scalar expectation values do not break symmetries in the large N limit; instead, the 
eigenvalues of the scalar operators are distributed over the transverse directions in a manner 
that respects the unbroken R-symmetries. The reasoning behind the boson rule is that a 
scalar expectation value allows an insertion of the operator dual to a supergravity fermion, 
generically of the form Yip, to deposit all of its momentum into the fermionic component 
ip, while the scalar Y is absorbed by the non- symmetry-breaking condensate. The results 
of m are clearest in cases where at least one of the independent chemical potentials is 


®Such a picture differs from the “gaugino” Fermi surfaces discussed in [7], as the Fermi surfaces in this 
case would be constructed from gauge invariant composite fields. Reconciling this interpretation with the 
scaling of the correlator and Luttinger’s theorem remains an interesting and unresolved issue. 
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absent in the black hole backgronnd. That is because unequal chemical potentials demand 
non-zero prohles for supergravity scalars whose held theory duals are expectation values 
of operators composed entirely from the held theory scalars whose non-symmetry-breaking 
condensates drive the reasoning behind the boson rule. It is unobvious how to extend the 
reasoning to the present case, where all four chemical potentials are equal, because then the 
non-symmetry-breaking supergravity scalars are altogether absent. It would be useful to 
examine supergravity constructions in which one, two, or three of the chemical potentials 
are turned oh, in order to try to ascertain whether some version of the boson rule can be 
applied even in the presence of a symmetry-breaking scalar like A. In the present case, it 


is interesting and suggestive to note from hgure 10 and equations (60)-(63) that the held 
theory operators which contribute to the fermion zero-energy mode at positive k are the ones 
whose bosonic components involve Y\ not Y^. For the Massive Boson domain wall, it would 
be in the spirit of the boson rule to speculate that this is because the deforming operator 
Os involves Y^ but not Y\. For the Massive Fermion domain wall, where the deforming 
operator Op involves only ijjA but not it is not clear how an argument in the style of 
the boson rule should go. We hope to report further on held theory interpretation in future 
work. 
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